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The  development  of  computational  methods  for  describing  the  adhe^vPpro^ss  was 
1  explored.  Serious  problems  currently  face  the  adhesive  technologist  in  that  some  of  the  best 
models  for  understanding  adhesion,  i.e.  Lewis  acid/base  neutralization  reactions,  are  based  on 
empirical  methods.  It  was  thought  that  computational  methods  may  show  a  way  for 
addressing  this  and  provide  a  first  principle  basis  for  understanding  adhesion. 

Current  limitations  in  software,  and  computer  hardware  do  not  yet  provide  for  using 
computers  as  a  common  tool  for  understanding  the  adhesive  process.  HOMO  and  LUMO 
values  are  thought  to  be  directly  identified  with  Lewis  acid/base  reaction  pairs.  TTicse  values 
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limited  value.  Further  work  in  our  understanding  of  quantum  physics  and  software  to  deal 
with  large  arrays  of  atoms  in  a  rigorous  manner  will  needed  before  computational  methods 
can  directly  speak  to  the  problems  of  adhesion  science. 
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Introduction 


The  goal  of  this  research  program  was  to  develop  a  set  of  computer- 
based  tools  to  help  understanding  the  adhesive  process.  To 
accomplish  these  tasks,  we  encountered  any  number  of  difficulties. 
These  difficulties  were  of  a  hardware/software,  personnel,  and 
scientific  basis.  While  we  learned  a  lot  about  the  nature  of  dealing 
with  the  subject  of  computer-based  research  in  adhesion  and 
adhesives,  circumstances  worked  against  our  meeting  the  specifics 
of  the  originally  defined  research  proposal. 


TASK  ONE:  Modify  the  already  existing  MACROMODEL  program  to 
make  it  more  suitable  for  adhesion  research. 


This  task  was  undertaken  with  the  understanding  that  we  would  have 
the  ongoing  help  from  the  director  of  the  Macromodel  program 
director.  Dr.  Clark  Still  of  Columbia.  Dr.  Still  had  done  all  past  work 
within  the  VAX  environment,  typically  using  Evans  and  Sutherland 
high  resolution  monitors,  equipment  we  have  available  to  us. 

Shortly  after  the  project  started.  Dr.  Still's  Laboratory  acquired  a 
Silicon  Graphics  platform.  All  work  on  the  further  development  or 
support  of  the  VAX-based  Macromodel  stopped  with  the  switch  to 
the  UNIX  environment.  Dr.  Still  informed  us  that  he  was  no  longer 
interested  in  any  aspect  of  VAX-based  Macromodel  development.  ,  , 

Since  we  did  not  at  that  time  have  access  to  an  SGI  on  campus,  we  - - - - - 

were  out  of  luck.  I  *y~ - 
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Another  aspect  of  this  task  was  the  hiring  of  a  post  doctorate  to  do 
the  programming.  I  quickly  learned  that  post-doc’s  in  computational 
chemistry  are  non-existent.  I  made  numerous  attempts  to  find  such 
a  person,  but  to  no  avail.  I  did  finally  find  a  person  with  a  dual 
bachelors  degree  in  computer  science  and  biology.  This  person,  a 
Chinese  student,  had  applied  for  a  position  with  the  VADMS 
(Visualization  and  Design  in  Molecular  Systems)  Center  here  at  WSU 
at  a  technical  programmer.  He  interviewed  second  out  of  9  and  was 
beat  out  by  someone  more  familiarity  with  genetics.  I  felt  fortunate 
to  have  such  a  person,  but  was  grossly  mislead.  This  person  was 
more  interested  in  getting  a  green  card  and  getting  married  to  a 
national  than  in  the  program  of  research.  I  had  to  release  him  (first 
time  in  my  professional  career  I’ve  ever  had  to  release  anyone.) 

Without  the  help  we  had  expected  form  Clark  Still’s  lab  and  without 
a  full  time  programmer,  we  were  not  able  to  meet  this  section  of 
the  program. 

Task  TWO.  Use  currently  existing  computational  tools  to  develop  a 
library  of  surfaces  for  general  use.  Compare  existing  computational 
derived  energies  of  adhesion  with  literature  values. 

I  had  some  serious  problem  with  this  topic.  First,  the  concept  of 
comparing  computational  values  with  values  derived  from  the 
literature.  Energies  are  calculated  in  typical  molecular  mechanics 
force  fields  as  kJ/mol  or  kCal/mol.  One  of  the  advantages  of 
Macromodel  as  a  research  tool  is  that  it  provides  within  the  shell, 
three  different  commonly  used  force  fields,  MM2,  AMBER,  and 
OPLS/A.  Each  force  field  has  a  unique  method  for  calculating  the 
internal  energy  of  a  system.  The  energy  from  one  set  of  calculations 
bears  little  resemblance  to  the  ehergy  calculated  with  another.  For 
instance,  propyl  alcohol  when  minimized  with  MM2  will  have  an 
internal  energy  of  10.22  kJ/mol,  AMBER  yields  3.01  kJ/mol,  and 
OPLS/A  yields  -0.01  kJ/mol.  This  is  a  seriously  wide  rahge  of 
values,  and  all  derived  from  within  Macromodel.  To  then  compare 
the  energies  derived  from  a  Macromodel  environment  with  another 
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modeling  environment,  such  as  BI06RAF,  or  SYBYL  quickly  became  a 
questionable  pursuit.  There  is  no  basis  for  holding  one  set  of  values 
in  higher  regard  than  another. 

A  similar  problem  existed  on  a  more  nuts-and-bolts  basis  when 
trying  to  develop  a  library  of  surfaces  to  use  for  modeling.  Given 
the  problems  of  different  energies  from  force  fields,  the  mere 
compatibility  of  one  program’s  data  files  with  another  became  a 
serious  point  of  consideration. 

Well  into  the  second  year  of  this  project,  Boeing  Aircraft  Comp, 
made  it  possible  for  us  to  purchase  used  SGI  monitors  for  a  nominal 
amount.  The  VADMS  laboratory  purchased  several  of  these  machines 
and  several  companies  were  invited  in  to  our  laboratory  to 
demonstrate  their  software.  While  some  software  claimed 
compatibility  with  Macromodel,  none  could  actually  demonstrate  it. 
Further,  the  major  commercial  programs  have  a  difficult  time  in 
reading  each  others  data  files,  not  surprising  considering  the 
competitive  nature  of  this  industry.  Everyone  was  “working”  on 
ways  of  providing  such  a  compatibility  in  the  future. 

TASK  THREE:  Development  of  a  force  field  for  computationally 
deriving  adhesion  forces. 

While  this  task  was  not  originally  expected  to  take  up  a  significant 
portion  of  the  overall  projects  time,  it  turned  out  to  be  one  of  the 
most  fruitful  and  one  where  significant  advances  in  understanding 
the  adhesive  process  took  place. 

There  is  an  understanding  that  Lewis  acid/base  or  donor/acceptor 
interactions  are  somehow  very  closely  identified  with  the  adhesive 
process.  Lee  ( 1 , 2)  in  two  book  has  developed  this  school  of  thought. 
The  work  of  Fowkes  (3,  4)  points  to  a  direct  relationship.  In  our  own 
work  we  solved  a  long-standing  mystery  of  wood-solvent  using  the 
donor  number  (DN)  and  acceptor  numbers  (AN)  of  Gutmann  (5).  Wood 
was  selected  because  it  is  a  natural  composite,  i.e.  cellulose  bonded 
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together  with  a  polyphenolic  material  lignin.  We  were  able  to 
explain  the  interactions  of  solvents  on  wood,  correlating  the  MOE, 
rupture  strength  and  acoustic  emission  of  wood  under  bending  loads, 
often  with  revalues  greater  than  0.99  with  simple  models  of  donor 
and  acceptor  number  (6).  A  particularly  strong  model  correlated 
strength  values  of  wood  against  DN/(mv)'^^  where  mv  is  the 
molecular  volume  of  the  treating  solvent.  Since  DN  are  measured  in 
calories  or  joules,  we  have  energy  divided  by  the  cube  root  of 
volume,  which  in  the  cgs  system,  is  cm.  Thus,  energy  divided  by  cm 
is  force.  We  have  determined  a  very  strong  correlation  between  the 
physical  reaction  of  a  polymeric  solid  and  a  force-related  value  of 
the  treating  solvent. 

The  general  problem  with  this  approach  is  that  the  acid/base 
numbers  or  donor/acceptor  numbers  are  all  derived  using  empirical 
methods.  There  is  no  first  principle  way  of  calculating  donor 
numbers  and  accefftor  numbers.  For  adhesive  researchers  this  is 
most  frustrating. 

We  spent  a  good  deal  of  time  focusing  in  on  this  problem.  Our  initial 
approach  was  to  calculate  the  highest  occupied  (HOMO)  and  lowest 
unoccupied  (LUMO)  molecular  orbitals  for  molecules  and  then 
attempt  to  determine  if  we  could  identify  relationships  between 
these  values  and  the  donor  number  and  acceptor  numbers  of  Gutmann. 

Two  approaches  were  considered.  First  we  were  able  to  get  a  little 
time  on  the  super  computer  platform  at  Cornell  and  ran  a  series  of 
ab  initio  calculations  using  a  variety  of  basis  sets  for  a  variety  of 
solvents.  On  local  computers  we  did  similar  runs  using  MOPAC,  a 
semi-empirical  software  program. 

Table  1  shows  the  results  of  a  series  of  runs  made  with  Gaussian  86. 
Methanol  was  selected  as  a  target  molecule  because  it  is  a  well 
studies  solvent  in  Lewis  acid-base  reactions.  Here  we  see  that 
selection  of  various  basis  sets  generated  energy  gaps  ranging  from 
16.340  to  25.368  ev  or  a  difference  of  55.25  %.  Table  2  shows  a  set 
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of  data  for  methanol  determined  with  MOPAC.  Here  we  see  a  closer 
range  of  energy  gaps,  but  all  values  falling  below  the  lowest  for 
methanol  determined  with  Gaussian  86. 

A  series  of  docking  experiments  was  conducted  using  MOPAC. 

Docking  is  defined  as  the  minimization  of  two  or  more  molecules  in 
close  proximity,  mimicking  the  adhesive  process.  Either  methanol  or 
urea  was  docked  against  graphite  modeled  with  a  varying  number  of 
unsaturated  rings.  Table  3  shows  a  range  of  68%  for  energy  gaps  of 
methanol  against  graphite  while  Table  4  shows  a  range  of  values  of 
277%  for  urea  against  graphite.  The  energy  gaps,  related  to  HOMO 
and  LUMO  values  and  thus  theoretically  related  to  acid/base  values 
is  very  sensitive  to  the  size  of  the  graphite  surface  used.  To  explore 
this  phenomenon,  a  series  of  runs  were  made  varying  the  size  of  the 
graphite  substrate  with  MOPAC.  Table  5  shows  that  energy  gap 
values  vary  by  as  much  as  9 1  %  for  the  larger  surfaces.  Since  we 
were  limited  by  the  size  of  our  computer  is  was  not  possible  to 
explore  significantly  larger  systems. 

The  outcome  from  these  computer  experiments  was  a  bit  frustrating 
from  two  points  of  view.  First,  we  could  identify  any  number  of 
legitimate  parameters  that  would  yield  most  any  HOMO  and  LUMO 
values  desired,  i.e.  determine  the  HOMO  and  LUMO  values  desired  for 
any  molecule,  we  can  tell  you  how  to  “legitimately”  calculate  them. 
This  is  a  problem  with  our  current  understanding  of  quantum  physics 
and  related  limitations  with  software. 

Another  aspect  of  this  problem  was  in  any  attempt  to  correlate 
these  values  with  calculations  based  on  molecular  mechanics  force 
fields.  For  example,  the  reasonable  approach  was  to  create  a  model 
graphite  surface  and  then  dock  a  molecule  such  as  methanol  or  urea 
to  it.  By  determining  the  total  energy,  and  then  the  energy  of  each 
component,  it  should  be  possible  to  determine  the  secohdary  bond 
energies  which  are  unique  to  the  total  assembly  which  are  not 
present  with  individual  components.  Such  an  experiment  using 
molecular  mechanics  force  field  MM2  yield  secondary  forces 


5 


between  graphite  and  methanol  of  -0.0149  ev  and  graphite  and  urea 
of  -0.243  ev.  These  compare  within  an  order  of  magnitude  for  some 
values  of  the  energy  gap  for  similar  systems  when  determined  by 
MOPAC,  but  not  for  others.  Similarly,  MOPAC  will  agree  with  some 
values  for  the  energy  gaps  for  graphite  with  Gaussian  86  but  not 
others.  It  Is  heartening  to  note  that  these  values,  roughly  bracket 
the  range  of  energies  thought  to  be  responsible  for  adhesive 
systems,  1  ev  =  23.06  kCal/mole.  The  problem  remains.  There  Is  no 
obvious  way  to  determine  which  values  more  closely  represent 
reality.  We  are  limited  by  the  size  of  the  arrays  we  can  compute  and 
are  unable  to  work  with  realistically  large  systems. 

The  second  problem  Is  more  fundamental  In  our  understanding  of  the 
science  Involved.  Donor  numbers  as  published  by  Gutmann  (3)  can  be 
shown  to  directly  correlate  with  the  properties  of  solid  materials 
and  with  the  joint  strengths  of  adhesive  substrate  systems.  Donor 
numbers  are  determined  by  measuring  the  heat  of  Interaction  of  a 
solvent  against  a  known  acid,  typically  antimony  pentachlorlde.  We 
perceive  this  to  be  a  reaction  of  SbCls  acting  like  the  perfect 
acceptor  and  the  unknown  donating  In  some  measurable  fashion,  a 
share  In  an  electron  pair.  Acceptor  numbers  are  based  on  shifts 
within  the  ^'P  NMR  spectra  for  triethylphosphlne  oxide,  a  standard 
reference  base,  dissolved  In  various  solvent.  This  Is  how  we 
perceive  the  reaction.  In  fact  that  is  not  the  only  thing  that  Is 
happening  to  generate  heat  or  induce  the  NMR  shift.  In  the 
calorimeter  or  NMR  we  have,  in  addition  to  this  simplistic  reaction, 
a  contribution  of  electrostatic  forces,  London  forces,  hydrogen 
bonding,  and  other  forces  all  Influencing  the  reaction.  Thus: 

AHtotal  -  AHactd-base  AHLondon  AHeiectrostatic  ■’’AHhydrogen  bonding  AHoiher 

To  date  researchers  have  been  trying  to  computationally  derive  the 
energy  of  acid-base  interactions,  qz  dispersive  forces 
electrostatic  Interactions  and  call  that  value  the  forces  of  adhesion. 
Any  real  system,  something  as  simple  as  placing  a  drop  of  water  on 
a  clean  glass  slide  will  In  all  likelihood  Involve  acid-base 


neutralization  factors,  along  with  dispersive  interactions,  possibly 
some  chemical  bonding,  hydrogen  bonding,  and  perhaps  even  other 
factors  we  don’t  yet  know  of.  Alas,  this  project  has  demonstrated 
to  me  the  immense  complexity  of  defining  the  forces  of  adhesion 
computationally. 
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Table  I  GAUSSIAN  calculations  of  methanol  (orbital  energies) 


Basis  sets 


Energy 

level 

ST0-3G 

3-2 IG 

4-316 

6-3  IG 

6-31  IG 

1 

-20.2633 

-20.4353 

-20.5103 

-20.5504 

-20.5437 

2 

-1  1.0858 

-1  1.21  17 

-1  1  2523 

-1  1.2738 

-1  1.2750 

3 

-1.2785 

-1.3346 

-1.3515 

-1.3544 

-1.3576 

4 

-0.8770 

-0.9271 

-0.9372 

-0.9389 

-0.9396 

5 

-0.61  10 

-0.6692 

-0.6828 

-0.6840 

-0.6866 

6 

-0.5713 

-0.6139 

-0.6217 

-0.6218 

-0.6230 

7 

-0.5172 

-0.5771 

-0.5889 

-0.5897 

-0.5918 

8 

-0.4187 

-0.4809 

-0.4878 

-0.4888 

-0.4917 

9 

-0.3594 

-0.4345 

-0.4446 

-0.4456 

-0.4479 

10 

0.5729 

0.2661 

0.2274 

0.2224 

0.1526 

1  1 

0.6507 

0.3158 

0.2744 

0.2724 

0. 1844 

12 

0.6969 

0.3452 

0.3187 

0.3159 

0  2232 

13 

0.7234 

0.3554 

0.3233 

0.3225 

0.2303 

14 

0.7345 

0.4245 

0.3914 

0.3740 

0.3325 

15 

0.9419 

0.8503 

0.7662 

0.4567 

16 

0.9623 

0.8681 

0.7828 

0.4812 

17 

0.9677 

0.8723 

0.8014 

0.5980 

18 

1.2479 

1.1415 

1.1247 

0.6706 

19 

1.3147 

1.1860 

1.1735 

0.7316 

E-gap  (au) 

0.9323 

0.7006 

0.6720 

0.6680 

0.6005 

E-qap  (ev) 

25.368 

19.063 

18.285 

18.176 

16.340 

Table  2.  MOPAC  calculations  of  methanol:  orbital  energies  (ev) 


Energy 

Levels 

MINDO/3 

Hamiltonians 

MNDO 

AMI 

PM3 

1 

-32.0896 

-41.9264 

-37.8433 

-38.1971 

2 

-24.6268 

-27.8236 

-26.9845 

-26.4865 

3 

-17.3034 

-18.6721 

-18.2860 

-18.0183 

4 

-15.2853 

-15.3874 

-15.4416 

-15.5270 

5 

-14.0382 

-15.3171 

-15.2892 

-15.5051 

6 

-1  1.9328 

-12.8217 

-12.5023 

-12.4843 

7 

-1  1.0781 

-11.4124 

-1  1.1306 

-1  1.1294 

8 

2.1020 

3.7950 

3.7758 

3.5152 

9 

3.4907 

3.9294 

4.0984 

3.9023 

10 

3.7895 

4.6937 

4.2238 

4.21  1  1 

1  I 

4.0557 

5.0164 

4.5217 

4.3928 

12 

5.6273 

6.7591 

6.2475 

5.6422 

E-gap  (ev) 

13.180 

15.207 

14.906 

14.645 

Table  3.  MOPAC  calculations  of  Graphite  docking  with  methanol 
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